Nonlinear I-V, High Flux-Flow Velocity Instability and Evidence for Quantum 
Fluctuations in BSCCO Superconducting Films 
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The large flux-flow velocity behavior of the high temperature superconductor BSCCO in the mixed 
state is studied by measurement of the current-voltage characteristic as a function of temperature 
and magnetic field direction. We find strong evidence for vortex-antivortex pairs that are excited 
by quantum fluctuations and contribute to the dissipation when dissociated by a large current. 
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Type II superconductors in the mixed state exhibit 
non-linear current-voltage dependence associated with 
the long energy-relaxation time of the conduction elec- 
trons At high flux-flow velocities this non- linear be- 
havior gives rise to an instability which manifests itself 
as a distinct break in the current-voltage characteristic. 
This behavior originates from the nonequilibrium distri- 
bution of quasiparticles in a superconductor driven by 
an electric field [El and was first predicted by Larkin and 
Ovchinikov (LO) In the LO theory, as the vor- 

tex velocity increases, quasiparticles exit the core 

1 and 

the viscous drag is reduced, driving more quasiparticles 
from the shrinking core. When all quasiparticles leave 
the vortex core the resistance of the sample jumps to a 
value close to the normal state resistance. The voltage, 
V*, at which this instability occurs defines a critical ve- 
locity for the vortices, v* , independent of the magnetic 
field, H, via the relation: 



V* = -(v* x H) L 



(1) 



where L is a vector along the direction of the current 
flow whose magnitude is equal to the length of the cur- 
rent path. This phenomenon has been observed in both 
low-T c and high-T c superconductors. Since 

the main assumption in the LO theory is that the distri- 
bution of quasiparticles is spatially uniform, the theory 
may be valid only near T c or at relatively high fields 
where the field is homogeneous in the sample. This re- 
quires that the inelastic scattering length of the quasi- 
particles at the critical velocity be larger than the inter- 
vortex distance, l H @: (O^) 1 / 2 > l H = i%/H) 1/2 ] 
here D is the quasiparticle diffusion constant, Ti n is the 
inelastic scattering time, and <I>o is the flux quantum. In- 
deed, Doettinger et al. [[Io) estimated that for the case 
(Dt^) 1 / 2 <C Ih = (*o/i0 1/2 , the result for the critical 
velocity will be modified, yielding a magnetic field depen- 
dence for which u* =| v* |~ ($ / 'H) 1 / 2 /r in . While this 
formula seems to fit the data in the range of fields studied 
in (!(]], new lower field data shows strong deviations. 

In the present paper we investigate further the I — V 
characteristics of BSCCO single crystal films in the tem- 
perature range T c /2. By using both field dependence and 



angular dependence at a fixed magnetic field we confirm 
the vector product nature of equation [j]. Measuring at 
very low perpendicular fields, we show that V* across 
the sample is finite and independent of magnetic field 
down to H = with very weak temperature dependence. 
This observation indicates the existence of a large num- 
ber of quantum mechanically excited vortex-antivortex 
pairs, detected in our experiment by current induced dis- 
sociation of the pairs followed by flux flow. 

The samples were 300 A thick BSCCO films grown on 
(100) LaA103 using layer-by-layer molecular beam epi- 
taxy, with resistively measured T c 's of 82 K and a resis- 
tance linear in T above T c . The films were c-axis ori- 
ented, hence I — V measurements represent properties 
of the a-b plane. Sample composition was carefully con- 
trolled to ensure the absence of second-phase defects Q . 
High-resolution photoemission on films grown under the 
same conditions show normal and superconducting spec- 
tra consistent with the bulk-crystal analog |l2| . 

A split-coil magnet was used to allow the variation of 
the angle between the plane of the sample and the field 
with a resolution of 0.05 degrees. Experimental proce- 
dure consisted of single triangular wave sweeps of the 
current and simultaneous measurement of the voltage 
across the sample. Sweep rates were varied over ^1.5 
orders of magnitude with no apparent effect to the sam- 
ples' behavior; this is evidence that sample heating did 
not occur at these sweep rates. Also, the self-field at the 
highest currents used is too small to explain any of the 
results discussed below. Figure 1 depicts a typical sweep 
at T = 35 K where the instability is clearly visible and 
strong hysteresis is observed. In what follows the values 
of V* and /*, the instantaneous current at which the in- 
stability occurs, were measured on the rising side of the 
current pulse. Also apparent in figure 1 is the strong 
non-linear regime before the instability. This nonlinear 
behavior will be discussed later. 

As noted above, the instability voltage was measured 
for several values of the magnetic field varying the an- 
gle 9 between the plane of the sample and the magnetic 
field from to 90 degrees at various temperatures. Fig- 
ure 2 depicts a set of curves taken at three different field 
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magnitudes as a function of the angle. Multiplying the 
applied field by sinO, we plot in the inset of figure 2 
the instability voltage as a function of this perpendicular 
component only. The collapse of all points onto a single 
curve in the whole range of fields up to 2.0 T is direct 
evidence for the validity of the triple-product in equation 
|l|. It also suggests that V* is controlled by vortex dissi- 
pation. Through variation of the angle of the field, very 
small values of the perpendicular component of the field 
can be applied to the sample. As the point where V* is 
minimum is approached, a remarkable result is obtained: 
V* is not zero at H± = Hsin9—Q. In fact, careful exam- 
ination of the results presented in ref. jl0| suggests that 
this phenomenon is also observed in YBCO thin films at 
sufficiently low temperatures (this is manifested by the 
fact that plotting u* as a function of H^ 1 ' 2 shows devi- 
ations towards saturation at low fields). 

The observation of a finite V* at low fields suggests 
that the extraction of a critical velocity from the V* data 
by using equation [l] (i.e. dividing the V* data by the 
perpendicular component of the field) may obscure the 
important features at low fields. However, it is evident 
from figure 2 that past a well-defined low field, V* indeed 
increases sub-linearly as a function of the field's perpen- 
dicular component similar to what is suggested by Jl0| . 
We find that the most striking way to present the data 
is by plotting V* as a function of the square root of the 
perpendicular component. This is done in figure 3 where 
we present data for three temperatures in the range of 
~ T c /2. The data clearly show a linear dependence of 

1 /2 

V* in H/ at higher fields that saturates to a finite V* 

1/2 

value at low fields. It is interesting to note that the H± 
dependence extrapolates to V*=0 from this higher field 
regime. In other words, at higher fields the system does 
not know that it will have a finite value of V* at low fields. 
Also, in no region of fields was there found a crossover 
to V* linear with field as predicted by the LO theory 
0. While for lower temperatures the saturation value 
decreases, it is evident from the data that this saturation 
remains finite even at T = K. 

The above results suggest that vortex dissipation is 
observed at low temperatures and zero magnetic field. 
Such an effect can be due to vortex-antivortex pairs that 
are excited by some mechanism. Because the thermal 
energy is low at ~ Tkt/^, thermally excited vortex pairs 
are excluded p3| . At 25 K e.g. we expect a density less 
than 10 -7 of thermally excited vortices. Moreover, the 
sample is very clean and pairs created by static disorder 
can be ruled out. This can also be determined from the 
data obtained at higher fields: V* extrapolates to zero 

1/2 

(no dissipation) as H± tends to zero. For a dirty sample 
with pairs created by disorder it is expected that this 
intercept remains finite. Thus, we conclude that these 
are observations of dissipation due to vortex-antivortex 
pairs that are excited by quantum phase fluctuations. 
These are quenched at higher fields, thus recovering the 
LO regime. The mechanism that causes dissipation is 



therefore the dissociation of these pairs. 

In order to understand the low current data we first fit 
the observed voltage as a power of the current 



v(i) = v [(i/i c r i] 



(2) 



The nonlinear regime before the instability is not hys- 
teretic and fits well a power law above I c giving a = 
z + 1 = 2.5 (while a — 2.5 gives best fit, a reasonable 
fit can be obtained with a in the range 2.4 to 3.0), and 
I c = J c wd ~ 0.008 A is a critical current associated with 
the onset of strong nonlinerarity due to vortex pair for- 
mation (w is the width of the sample, d its thickness). 
Here we introduce the dynamical exponent z, which is 
found to be smaller than 2. Note that a is equal to 3 
(i.e. z — 2) in classical Kosterlitz-Thouless theory pj 
close to Tkt 17 1- Using the data presented in figure 1, 
figure 4 depicts V vs. I. This figure shows that there 
is an appreciable deviation from equation || for / near 
I c indicating the presence of a population of pairs with 
spacing greater than the average spacing. To understand 
I c and its relation to the average spacing between vor- 
tex pairs, we note that the energy of a vortex pair with 
separation r in the presence of a current density J is: 



E(r) = U v ln(r/£) + 2E C - (l/c)Jd^ r 



(3) 



where U v — AksTxr [|15[ is the Kosterlitz-Thouless tran- 
sition scale and E c = 0.39U V |Q is the core energy. 
Minimization gives r av = cU v / J c d$Q. Thus, I c found 
in figure 4 implies r av ~ 105 A, or a size equivalent 
to a few vortex cores. If we were observing a thermal- 
driven Kosterlitz-Thouless transition p6| , the non-linear 
I — V behavior before the instability would originate 
in thermally-activated hopping over the barrier E(r av ). 
Choosing £ = 20 A in equation ||, this rate can be calcu- 
lated for the temperatures at which the data was taken. 
We find that at 55 K the rate is suppressed by a fac- 
tor 5xl0 -4 , while at 25 K it is suppressed by a factor 
6xl0 -8 . Moreover, we note that between 25 K and 35 
K the rate would change by a factor of 100, implying a 
similar change in the leveling field, while experimentally 
this change is less than a factor of two. We therefore con- 
clude that the observed nonlinear behavior and voltage 
instability must originate from vortex-antivortex pairs in- 
duced by quantum fluctuations and not from thermally 
created pairs. The smaller dynamical exponent can also 
point to a quantum mechanical dominated result as all 
classical results on high-T c show increased z above 2 [jlB] . 

To model the distribution of pair spacings we assume 
a simple harmonic oscillator model for the wave function 
of the vortex-antivortex pairs: 



ip(r) 



-i(r/a) 2 



(4) 



where a ~ r av is of order a few vortex core distances. 
Then we can relate the distribution of pair spacings to 
the distribution of critical currents by writing for fixed r 



2 



v(r,J) = p[(J/J c (r)r-l]6(J-J c (r)) 



(5) 



where p is the density of pairs. Using equation ^| the 
voltage response as a function of r can then be obtained 
by integrating over the pair distribution as 

/>oo 

V(J) (xp dr\*jj(r)\ 2 v(r,J) (6) 

J A/ J 

which may be written in dimcnsionless form as 

/>oo 

V(J) =V [{Jx) a - l)]x 2 e- x2 dx (7) 

J A/ J 

Here J = J/J c (a), A = (cwU v / irha) , and eVo = 
(l/c)pJd<& a represents the voltage scale corresponding 
to dissipation at a length scale a. The resulting J depen- 
dence is parameter free and depends only on the scal- 
ing index a. The result for a — 2.5 is shown in fig- 
ure 4 as an inset and fits the data remarkably well. We 
take this as an indication that the low current tail in the 
I — V curves is a direct measure of the distribution of 
vortex pairs excited by zero point motion. Coupled with 
the high current instability showing that the resistivity 
is due to vortex motion, we believe this is the first di- 
rect evidence of quantum excitation of vortex pairs in 
high-T c superconductors at temperatures well below T c . 
This lends support to the idea that the loss of phase co- 
herence in high-T c superconductors may be directly at- 
tributed to quantum-fluctuation-induced winding of the 
superconducting order parameter |l9|,^o| . 

Observation of large dissipation in high-Tc films at low 
temperatures was previously reported by Corson et al. 
pl| , who measured complex conductivities in the fre- 
quency range of 100-600 GHz of BSCCO films from the 
same source as those used in this work. At low temper- 
atures they found a residual real-part conductivity that 
scales with the T = superfluid density and is larger 
than any contribution that could be due to thermal pop- 
ulation of quasiparticles at the d-wave gap nodes. 

To ascertain the relevance of the superconductor pair- 
ing symmetry to our conclusions, we performed similar 
measurements on 120 A film of amorphous MoGe previ- 
ously used to study vortex dynamics by Hellerqvist et al. 
p2j . These films were chosen for their very short coher- 
ence length, large penetration depth, and relatively weak 
pinning, similar to very thin high-T c films . The ob- 
served behavior of the MoGe films resembles that of the 
BSCCO films; the inset in figure 3 shows a V* curve at 
T = 4.3 K (T c = 6.6 K). Again, the extrapolation of V* 

as a function of H± 2 gives zero dissipation at zero field. 
However, it is clear that this system is not a d-wave super- 
conductor, thus the dissipation is wholly due to vortices. 
Note also that the density of vortex-antivortex pairs is 
rather high. The low temperature value can be extracted 
from the crossover field. At T = 25 K this field is ap- 
proximately 0.25 T, of order 0.5% of H C 2, indicating a 
pair density of similar value. The application of a finite 



field that creates vortices with only one sign therefore 
suppresses the quantum fluctuations as expected. 

Finally, we would like to comment on the consequences 
of the H 1 / 2 dependence of V* . It is interesting to con- 
nect our findings to other measurements related to the 
size of the vortex core, a key feature in the LO theory. 
Note that LO's result is derived from the calculation of 
the velocity at which the viscous- friction force has a max- 
imum. Their expression for the critical velocity has the 
form v* 2 ~ D/T in ~ £, 2 /(T m T GL ), where t gl is the fun- 
damental Ginzburg-landau time. As noted above, at low 
fields u* ~ {^o/B) 1 / 2 1 T in . Thus we may conclude that 
in this regime £ ~ Ih ~ H^ 1 / 2 and that this apparent 
increase in vortex core size is a consequence of the fact 
that the quasiparticle distribution shows a strong spa- 
tial inhomogeneity with portions of the superconducting 
phase still displaying the equilibrium distribution. Simi- 
lar behavior, an apparent increase in vortex core size, was 
also observed in recent muon-spin-rotation spectroscopy 
experiments on hig h- T c [ p4| and conventional layered su- 
perconductors |25[] , as well as in specific heat measure- 
ments |2f|. We believe that all these measurements are 
in fact observing the strongly inhomogcneous quasiparti- 
cle distribution in the system at low fields. 

In summary, we present evidence that in thin films of 
BSCCO, at low temperatures and low magnetic field (in- 
cluding H = T), there exists anomalous dissipation 
that is similar to vortex dissipation observed at finite 
magnetic fields. This provides direct evidence for the ex- 
istence of vortex-antivortex excitations due to quantum 
fluctuations of the phase of the superconductor. 
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FIG. 1. Typical I-V measurement cycle. The current 
is ramped in a triangular fashion and the corresponding 
voltage is recorded. V* marks the onset of the instability 
(see text). 

FIG. 2. Instability voltage as a function of angle for 
several magnetic fields at T = 35 K. Inset shows data 
collapse when plotted against the perpendicular (to the 
sample) component of the magnetic field. 

FIG. 3. Instability voltage vs. H 1 / 2 at two different 
temperatures (T c = 82 K). Inset shows similar data for a 
120 A sample of MoGe (see text). 

FIG. 4. The data of figure 1 plotted as V vs. /. Insert 
shows the fit to equation I?]. I c is discussed in the text. 
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